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Abstract 

Boundary correlation functions of the six and nineteen vertex models on an N x N lattice with 
domain wall boundary conditions are studied. The general expression of the boundary correlation 
functions is obtained for the six vertex model by use of the quantum inverse scattering method. 
For the nineteen vertex model, the boundary correlation functions are shown to be expressed in 
terms of those for the six vertex model. 
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1 Introduction 

The six vertex model is one of the most fundamental exactly solved models in statistical physics 
[H |2j [3] |4] . Not only the periodic boundary condition but also the domain wall boundary condition 
is an interesting boundary condition. For example, the partition function is deeply related to the 
norm [5] and the scalar product [6 of the XXZ chain. The determinant formula of the partition 
function [7J[5] lead Slavnov [5J to obtain a compact representation of the scalar product, which plays a 
fundamental role in calculating correlation functions of the XXZ chain [SJ HU H2] • The determinant 
formula also led to a deep advance in enumerative combinatorics (T3l [T4J [15] . For example, it was 
used to give a concise proof of the numbers of the alternating sign matrices for a given size. Recently, 
the correspondences between the partition function and the Schur polynomial [IB] and KP r function 
[17j have been revealed. The determinant representations of partition functions have been extended 
to other models such as the higher spin vertex models [IB], Felderhof models [TS] and so on. The 
domain wall boundary conditions are also interesting from the physical point of view since it exhibits 
phase separation phenomena [201 HU HH H31 HI] ■ 

The calculation of correlation functions are also interesting in the domain wall boundary condition 
itself. Several kinds of them such as the boundary one point functions, two point functions, boundary 
polarization [2B1 123 HE] and the emptiness formation probability [5S] have been calculated. 

In this paper, we calculate boundary correlation functions for the six and nineteen vertex models 
on an N x N lattice with domain wall boundary condition. For the six vertex model, we use the 
quantum inverse scattering method, and apply the approach of 29J to build and solve two recursive 
relations for the boundary correlation functions, providing for them a general expression. The bound- 
ary correlation functions we consider includes the boundary polarization and boundary emptiness 
formation probability (EFP) as special cases. 

Next, by use of fusion, we show that the boundary correlation functions for the nineteen vertex 
(Fateev-Zamolodchikov [3D]) model can be reduced to those for the six vertex model. In particular, 
the EFP of length s for the nineteen vertex model reduces to that of length 2s for the six vertex 
model. 



The outline of this paper is as follows. In the next section, we define the six vertex model with 
domain wall boundary condition. The general expression for the boundary correlation functions is 
obtained in section 3. In section 4, the boundary correlation functions of the nineteen vertex model 
are considered by use of of fusion. The emptiness formation probability for the nineteen vertex model 
is expressed in the determinant form in the homogeneous limit in section 5. 



2 Six vertex model 



The six vertex model is a model in statistical mechanics, whose local states are associated with edges 
of a square lattice, which can take two values. The Boltzmann weights are assigned to its vertices, 
and each weight is determined by the configuration around a vertex. What plays the fundamental 
role is the i?-matrix 
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(i) 



which satisfies the Yang-Baxter equation 

i?i 2 (A, v)R lz {\ [i)R 2 3(v, = R23(v, n)Ri 3 {X, /x)i?i 2 (A, v). 



(2) 



We consider the six vertex model on a N x N lattice depicted in Figure. The spins are aligned all up 
at the bottom and right boundaries, and all down at the top and left boundaries. At the intersection 
of the a-th row (from the bottom) and the fc-th column (from the left), we associate the statistical 
weight 



where 



£ a fc(A a , v k ) = sh(A Q - v k + rj/2)R ak (\ a - 77/2, v k ) 
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a(A, v) = sh(A - v + rj/2), b(\, v) = sh(A - v - rj/2), c = shr]. 



(3) 



(4) 



We refer to the a-th row as the auxiliary space V a and the fc-th column as the quantum space TL k . Let 
us denote {A} = {Ai, A 2 , . . . , Ajv}, {^} = {^1, ^2, • ■ • , ^n}, and the basis (dual basis) of the spin-1/2 
representation as |+), |— ) ((+|, (— |). 

The partition function of the six vertex model, which is the summation of products of statistical 
weights over all possible configurations can be formally represented as 
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where ||+) = ®f =1 |+) fe , ||-) = ®£ijj-}fc, <+|| = <8>£L lfe <+|, (-|| = ®f =lfe (-|, and we distinguish the 
spins on the quantum and auxiliary spaces by the subscripts " g" and " a" . The partition function has 
the following determinant form [8] 
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Figure 1: The six vertex model with domain wall boundary condition. 



where 



d(X,v) = sh(X-v), M ak = (p(X a ,Vk), 



a(X, v)b{X, v) 



Introducing the monodromy matrix 

T a (X a , {v}) =C a N\Xci,VN) ■ ■ ■ £al(X a > V\) 

_( A(X a ,{u}) B(X a ,{u}) 
\ C(\*,{v}) D(X a ,{u}) 

the partition function can be represented as 

Zn({X}, M) = q (-\\B(X N , M) ■ • ■ B(Ai, M)||+) a 

From the Yang-Baxter equation, one has 

R a p{p. - X)T a (ji, {iy})T p (X, {v}) = T p (X, {v})T a (n, {v})R al) {p, - A). 

From (ITU1) . one has 

A(X, {v})B{^ {v}) = /(A, n)B(ji, {v})A(X, {v}) + g(ji, X)B(X, {v})A{^ {v}), 
B(X, M)A(fi, {v}) = /(A, fJ,)A(p,, M)B(X, M) + g(fi, X)A(X, {v})B{^, {„}), 
B(X, M)B(/x, {v}) = B(p, M)B(X, {i/}), 



where 
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3 Boundary correlation functions of the six vertex model 

In this section, we consider the following boundary correlation functions 
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where irt = \+}kk(+\ and irl = \—}kk{~\ is a projection onto the up and down spin respectively. Some 
special cases of this general boundary correlation function reduces to the ones previously considered 
[2"51 [23 [571 E2j • We calculate the boundary correlation functions by use of the quantum inverse 
scattering method, applying the approach of [25] . 



A 



N 



X 



r+1 

Ai 



Figure 2: An example of boundary correlation function. 
First, note that (HI)]) can be expressed as 

j£' ei '"' ' £s) ({A}, M) = q {-\\B{\ N , M) ■ ■ ■ B(X r+u M) n 4^(A,, M) • • • B(Ai, M)||+)„ (17) 

fe=i 

in the quantum inverse scattering language. 

We introduce the following two-site model |10| in order to obtain recursive relations between 
boundary correlation functions of different lattice sizes. 

T(A, M) =T 2 (A, {^VOT^A, ^i) (18) 
T 2 {\, j>}Vi) ^C aN (X, v N ) ■ ■ ■ C a2 (X, v-x) 



G 2 {\{v}\vx) ^ 2 (A,MVx) y (19) 
T 1 (A,i/ 1 )=£ al (A,i/ 1 ). (20) 
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Applying 

i(+|B(A, M)|+>i = b(X, u 1 )B 2 (X, {u}\vi), 
1 (-\B(X,{v})\+) 1 = cA 2 (X,{v}\u 1 ), 
i(+| J B(A,H)|-)i = > 

i(-|B(A, W)h)i = o(A, vi)B 2 (X, MVi), (21) 

iteratively, one has 

71 

! (+|fl(A„, M) • • • 5(Ai, M)|+)i = J] 6(Aj, i/i)S 2 (A n) MW) • • • B 2 (X 1 , M W), (22) 

71 

i(-\B(X n , {u}) ■ ■ ■ B(X U M)|-)i = J] o(Ai,n)B a (A„, WW) • • • B 2 (Xi, M\n), (23) 
1 (-|S(A„,M)---B(A 1 ,M)|+) 1 

n n a—1 

= ^ a(A l s,J/i)c JJ 6(A l 8,^i)-B 2 (A„,{!/}\^i) • • ■ B 2 (X a+ i,{v}\vi) 

a=l/3=a+l /3 = 1 

x A 2 (X a , {i/}\^)S 2 (A a _i > {^W) ■ ' ' B2(Ai, MVi). (24) 

Combining 

A 2 (A, M\^)5 2 ( M , {^}\^) =f(X,fx)B 2 (fx, M\^)A 2 (A, 

+ ff 0*, A)5 2 (A, {v}\ Vl )A 2 {n, {v}\vx), (25) 



and 



we get 



iV 



A 2 (A, {v}\vi) ®f =2 |+) fe = [] «*) ®k=2 \+)k, (26) 



fc= 2 



i<-|B(A r> H).-.B(A 1> H)||+> 

r r r jy r 

= e c n 6 (^,^i) n /( a - ^ n °( a - ^ n b ^ mw ®*u i+>*- ( 27 ) 

a = l 0=1 /3=1 fc= 2 fc=l 

In the same way as (|2T|). we can also show the following relation 
(-||B(A^,M)..- J B(A P+1 ,{i/})|+)i 

N N N N N 

= e c n «( a /3^i) n f(^^ a )i[b(x a ,^)^ =2k (-\ n b^a^hw), (28) 

a=r+l 0=r+l (3=r+l fc= 2 fe=r+l 

utilizing ([24|) . 

£? 2 (A, {y}\u 1 )A 2 {^ =f(X, f i)A 2 (fi, {v}\vi)B 2 (\, 

+ <?(//, A)A 2 (A, {i/}\^)B 2 (|i, (29) 

and 

2V 

®£U*H^(A, (i/jW) - J] b(X,v k ) ®£L 2 fe (-|. (30) 

fc= 2 



5 



From (|23[) and (|27[) . one can derive one recursive relation for the boundary correlation functions 
between different lattice sizes [2"9"] 

N r r r N 

jrp-,e 2 ,- , £s ) ({A}; {v}) = -Q a ( Xf}fUl) c [] 6(A^) J] /(A Q , A^g) J] a(A a , u k ) 

/3=r+l a=l /3=1 /3=1 fc=2 

x ^ : l,e 2 ,..,e s)({AnAc[;MW) _ (31) 

We can obtain another recursive relation from (l2"2l and (l2"5)l 

r N N N N 



^ + ,e 2 ,-, £s ) ({A}iM)= -Q 6(A ^ iyi) £ c JJ a(A/3)!/i) J| /(A^,A a )n & (Aa, 



/3=1 a=r+l ,3=r+l j3=r+l k=2 

xTpr-^\{X}\X a ,MVi)- (32) 

Solving these two recursive relations ([3~T]) and (1521 , one obtains the general expression for the boundary 
correlation functions as 

-4 r ' ei '-' es) ({A},M) 

dct m({a } , M) n n5=l a(A „ ^ n jL r+1 b(\ ? ,vj) ai ^, r J^ ir • ' ^ ^ 

a 2 ^tti a s ^ai,-,a a _i 
x (_ 1 )E 1 < J<fc < s x(afc,«i) + Efc =1 (afc-l-K£fc + l)/2) + i:fc = i(efc + l)(A r -fc)/2 
s 

xn^(^) II ^ e *(^a,.^a*,^.^)detM({A}\{A ai ,-.-,A Q .},{i/}\{i^,-. •,!/,}), (33) 

j=l l<j<k<s 

where S^' r = {1, • • • , r}, S+* r = {r + 1, • • • , N}, e(A, v) = sh(A - v + rj), 

rr- n , II^ie(Aff,A) n^+i^,A) 

lh=i b (\v k ) 

_ U r p=id{.X^X)Yf p=r+1 e{X,X p ) 

H r (A) - — iv , (35) 

lL=i a ( A ^fc) 
. , a(A a ^)a(A afc ,^-) 

B ^(A Qi , X ak ,V h V k ) = -r- - r , (36) 

E (A Qj , A Qfc , Vj ,iy k ) = 1 — ^ , (37) 

(38) 

E+-{X aj ,X ak , Vj ,v k ) = ""'XTT^ ^ 

and x(/3,a) = 1 for /3 > a and otherwise. The proof is given in the Appendix. As a special 
case {tj = —,j — 1, • • • ,s), The boundary correlation function reduces to the emptiness formation 
probability [29^ (cf. JTD] ) , which gives the probability of finding a sequence of all spins down of length 
s from the left boundary. 
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4 Nineteen vertex model 



In this section, we consider the nineteen (spin-1 or Fateev-Zamolodchikov) vertex model. The nineteen 
vertex model can be constructed from the gauge transformed spin-1/2 R- matrix 



R+(A,v) = ^(A)0 fc H J R, fe (A^)0Ti(A)^ 1 H, 



(40) 



where 0(A) = diag(l,e A ) and the projection operator 
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(41) 



The basis (dual basis) of the spin-1 representation |1), |0), — 1) ((1|, (0|, (— 1|) is given in terms of basis 
of spin-1/2 representation as |1) = |+) <g> |+), |0) = (1 + e^)" 1 / 2 ^} <g) |-> + e~ r '\-) <8> |+)), | - 1) = 
|— ) <%> \— ). The gauge transformed spin-1 i?-matrix can be constructed as [3"T1 1521 153"] 



R .JK ( z > w ) =PlK-\ 1 1KPlJ,1J-\Rtj,1K{ Z + ™)R2J,2K-l( Z + V,™ + V) 
X R2J-l,2K( Z ' W ) R 2J-l,2K-l( Z ' W + V)P2J,2J-lP2K-l,2K- 

The symmetric spin-1 i?-matrix can be obtained from i?^(z,w) by gauging out factors as 
IPjk&w) = ^j 1 (z)^- K 1 (w)R 1 J +(z,w)^j(z)^ K (w), 



(42) 



(43) 



where &(z) = diag(l, e 2 , e 2z ). 

For the nineteen vertex model on a N x N lattice with domain wall boundary condition, all spins 
are aligned +1 at the bottom and right boundaries, and —1 at the top and left boundaries. At the 
intersection of the a-th row (from the bottom) and the fc-th column (from the left), the statistical 
weight L^z^Wk) = R l ak (z a - rj/2, w k ) is associated. We also set L^(z a , w k ) = R^(z a ~ ^/2, w k ), 

1/2 1/2 - ! - -4- 

Wk ( z a> w k) = Rak(z a ~ rj/2, w k ), Lj k (z a ,w k ) = Rak( z <x - W 2 : w k) for later convenience. 
We consider the boundary correlation functions for this nineteen vertex model 
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where < fc = \S k )kk(S k \,S k = 1,0, -1 and ||1) = ®f =1 |l) fe , ||-1> = ®£Li|-l)fc, (1|| = ®fc=ife(l|, <-l|| = 

We show that the above boundary correlation functions can be reduced to those for the six vertex 
model calculated in the previous section. Instead of directly dealing with (|4"4"| . we consider 
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We also define F^f 2 ^ r ' ei '"' ' Cs ' > and F^ 2+ ^ r ' ei ' "' ' e '^ as well, replacing L 1 by L 1 / 2 and L l l 2+ , respectively. 
From one can see 

^ 1 '- A) ({4 ) W) = 4 !rAr ' A) ({4,W) 1 (so) 

since 

^ + ({4,W) = e^--- 2 ^(--"/ 2 )^({z},W), (51) 
^+<rA,--A> (WjM) = e ?^ (52) 

Thus, we can consider F^^'^ 1 ' ' Ss ^ ({z}, {w}) instead, which is easier to handle than F$ r ' Su '" A) ({z}, {w}) 
itself. 

We reduce Fl^~^ r ' Sl ' ({2}, {w}) to the boundary correlation functions of the six vertex model 
by use of fusion. The following relations are used. 

P2K-1,2K R 4j i 2k( Z + V, w ) R ij,2K-l( z + V,™ + l)4j-l,2ir(^ W ) R 2,J -1.2K -1 ( z > w + V) P 2K-X,2K 
= P2K-l,2KRtj,2K( z + V, w ) R 2J,2K-l( z + + V) R tj -1,2k( Z : > W ) R 2J-1,2K-1 i z > W + V), ( 53 ) 

P2J,2J-1 R 2J,2K( Z + V, w) R 2J,2K-l( z + + V) R tj-l,2K:( Z > W ) R 2J-l,2K-l( Z ' w + V) P 2J,2J-1 

= P 2 J,2J-1 R 2J,2k( Z + V, w ) R 2.J,2K-l( Z + + V) R 2J-1,2k( Z ' W ) R 2J-l,2K-l( Z ' W + V) , (54) 

P 2 = P, P\ ± 1) = |±) ® |±), (±1\P = (±| ® (±|, (55) 

4 k P2k-l,2k=P2k-l,2k £ C t-^AT-^2k, (56) 

where C*j + = Cz\ = C\_ = C°_ + = 1 and otherwise. 
First, utilizing (1551) . ([53)1 and ([55)1 . one has 

AT iV s 



4 +(rA '- A) ({4,w)=a<iii ? <-iiin p ^-i,2 fc n 7 i Q (z a ,w)n^2 fe -i, 2fe ^p 2fc _ 1 , 2fc 

fe=l a=r+l fe=l 

r 

n f a(«a,M)||->„||+)„ (57) 



X 

a = l 

where ||+) = 8™ |+) fc , ||-> - S^Rfc, (+11 = ®^i*<+l, HI - ®^i*H 

iV 

f Q (z a , {w}) = JJ i2a 2 2fe( Z " + »7, Wfe)i2a 2 2fe-l( Z « + '?> + V) 



fe=X 

X 

fe=l 



II ^-1,2*: (^ W fc) i 2a 2 i,2fc-l( 2 a : , ™fc + »?), (58) 



T Q (z a , {w}) =P2a,2a-ir Q (z a , {w}). (59) 

Applying ([55]) . we have 

N N 

F^^ r ' Sl ' " ' Ss \{z}, {w}) = a (l|| g ( — 1|| JJ P2k-l,2k II 

fe=l a.— r+1 

s 

fe=l 62is-l |E2k=± 



n T «(^,W)||-)a||+> g - (60) 



X 

a=l 
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Using (f55|) and f[55j) , one gets 

s 

Fj r +Mi, - A) (w,M)= e n^t-^^+ii^-ii n r„(2r a ,M) 



A' 



eii'" ,£2s = ± fc=l a=r+l 
H <l-l^t II f a(«a, {te})||->a||+>, 



X 
fc=l 



,e2s=± fe=l 

where {2} = {z x ,Zi + ri,z 2 ,z 2 + rj, ■ ■ ■ ,z N ,z N + r]}, {w} = {wi + r),Wi,w 2 + 77, w 2 , ■ ■ ■ , w N + r),w N }. 
As the simplest case, one has [T5] 

Z 1 +({z},{w}) = Z 1 2 / N 2+ ({z},{w}). (62) 

Thus we have 



C ,fI, " A, (w.W)= E n^-.^r^^'^fi^i-i*}). (63) 

eir" ,«2s=± fc=l 

From one can see 



C'^"'''"'^^^})^^''''^^}.^}), (64) 

since 

^({*}J©})=e?^-*->E^-.+^V»( {g}|{fi)})l (65) 
^ +(2r ' E1 '- ' E2<) ({^}, {»}) = e 2E "- ™- 2E «- ^+*^f 2r ^ - {»}). (66) 

Combining (|50|) . and one finally has 

4 (rA, "" A) ({*}.{«})= E n^_ l£2fc F 2 1 / 2 ( 2 ----)({z},M), (67) 

,«2s = ± fc = l 

which means that the boundary correlation functions for the nineteen vertex model on an N x AT 
lattice with spectral parameters {z},{w} can be reduced to those for the six vertex model on a 
2N x 2N lattice with spectral parmeters {z},{iD} . Note that F^ l " 2r ' ei '"' ' e2s \{z}, {w}) is exactly 

■^ajv' 61 '"' ' e2 ^ {w}) in the previous section since the corresponding Boltzmann weights are different 
just by an overall factor, which do not affect correlation functions. 

5 Homogeneous limit of the emptiness formation probability 

Let us consider the homogeneous limit of the emptiness formation probability (EFP) for the nineteen 
vertex model. As a special case of (|67[) . one has 

F^ l - r \{z},{w}) = 4 2 ^- f,) ({z},{w}), (68) 

i.e., the EFP of length s for the nineteen vertex model with spectral parameters {z}, {w} reduces to 
the EFP of length 2s for the six vertex model with spectral parameters {z} = {zi, z\ + rj, z 2 ,z 2 + 
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7/, • • ■ , z N , z N + r]}, {w} = {wi +T), w\, u>2 +Tf, W2, • • • 7 wn+V, w n}- Let us set Zj as Zj = z + £j. One 
finds that T^ X ~ )2S) ({z},{w}) can be expressed in the determinant form as 

r(2r,(-) 2s )/, i - Xi X 2 X 3 



Fin {{z}Aw}) = - — r .. det^^ r , (69) 

2N det M({z},{w}) Ili< j<k< 2 s e(z + tJ,z + tk) £l =-= ea .=o' 



where 



llj=i []l^i a ( 2 /3 +?7,Wj)a 2 (2/3,Wi)a(^,^ +v)\[Y\p=r+i b {zf3 + il,w o )b 2 {z 01 w 3 )b{z^w 3 + r])] 

(70) 

^ = IIj 2 ll [11/3 = 1 f(fg. f + £ j) e ( Z /3 + gvf + ej)][II?=r+l rf Ol> Z + ^jKfg + gvf + ^ 



Xa = I1 [II a ( z + e 2j-i> w fc) ]~[ a(2 + e 2 j-i,Wfe+r?)] 
j=l k=j k=j+l 

s — 1 s s 

x II [ II a ( z + e 2j,w k ) H a(z + e 2 j,w k +ri)] 
j=i fe=j+i fc=j+i 

s fe-1 fe-1 
X II [ II & ( Z + £ 2fc-l, Wj) b(z + £2/0-1, Uij + Jj)] 
fe=2 j=l j=l 
s fe— 1 /c 

x II \n b ( z + t2k,w j )\[b{z + e 2k ,w. J +T 1 )l (72) 
fc=i j=i j=i 



and 'J' is a 27V x 27V matrix whose (j, fc)-th block matrix element is given by 

exp(O^-i) exp^c^J \ 

ex P ((^ + fjja^j ex P ((^ + »?)a £2 j y ' 



for j = 1, • • • ,N,k = 1, - ■ ■ ,s and 



Lp{zj,w k +rf) ip(zj,Wk) 
<P{zj + V, Wk + V) <p{zj+V,Wk) 



(73) 



(74) 



for j = I,-- - ,7V,fc = s + 1,-- - ,iV. 

Now let us take the homogeneous limit by putting £j,Wj,j = 1 , • • • , TV to zero in the order w\ — > 
0, . . . , U7]v ->■ 0, £i 0, . . . , £n -> 0. We have 



l(r,(-) s ) -,-(2r,(-) 2s ) 5^1 , , , Y 2 Y 3 

f n = F?n = det 1p= —. r 

det m rii< i </ £ <2 S sh fe - e k + V) 



(75) 

£l=---=£2s=0 
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where 



= ( _ 1) ,N-.(, + l)/ 2{n » =i(JV _ j) ! } 2 d 2,N-, - faQ) 

1 a"(z + 7 ? ,0)6( Ar - r ) ;s (z + ?7,0)a 2 '' s (z,0)6 2 ( Ar - r ) s (z,0)a'- s (z-r/,0)&( Ar - r ) s (z-77,0)' ^ j 

= S h r (- Cj - + 2 ?? )sh jv (- £j + v)sh N - r (-ej) 

2 }=i sh Ar (e, +z-r;/2)sh w (e J +z-3ry/2) ' 1 j 

y 3 = II shS ~ j+1 ( z + £ 2j-i + V/2)& s - j (z + e 2j _! - 77/2) 

s-1 

x JJ sh s ~ J ' (z + e 2j - + r?/2)sh s - j (z + e 2j - - r?/2) 

S 

x JJ sh* -1 ^ + e 2fe _i - 7 ? /2)sh fe - 1 (z + e 2k -i - 3ry/2) 
fe= 2 

s 

II shfe_1 ( z + ^ - »?/2)ah fc (* + e 2fc - 377/2), (78) 



x 
fe=i 



mis a 2iV x 2iV matrix whose (j, fc)-th block matrix element is given by 

di +k - 2 tp{z-r),Q) di+ k - 2 (p(z,0) \ 
di+ k - 2 <p{z,0) dl +k - 2 <p{z + ry, 0) J ' 



(79) 



for j, k = 1, • ■ • , N, and tp is a 2iV x 27V matrix whose (j, fc)-th block matrix element is given by 

gj+k-2 gj+k-2 \ 

MvdZ k ~i)di+ k - 2 eMvdZ)di+ k ~ 2 ) ' (80) 
for j = 1, • • • ,N,k = 1, - ■ ■ ,s and 

+ fe - s -^(z-77,0) 0*+ fc —V(*,O) \ 

for j = 1, • • • ,N,k = s + 1, - ■ ■ ,N. 

6 Conclusion 

In this paper, we considered correlation functions for the six and nineteen vertex models on an N x N 
lattice with domain wall boundary conditions. For the six vertex model, we derived the general expres- 
sion for the boundary correlation function by solving two recursive relations obtained by the quantum 
inverse scattering method. The result includes the boundary one point functions and emptiness for- 
mation probability as special cases. 

For the nineteen vertex (Fateev-Zamolodchikov) model, by use of fusion, we have shown that the 
boundary correlation functions reduce to those for the six vertex model. In particular, the emptiness 
formation probability of length s for the nineteen vertex model on an N x N lattice reduces to that 
of length 2s for the six vertex model on a 2N x 2N lattice with appropriate spectral parameters. 

The correlation functions "off the boundary" can in principle be expressed as a linear sum of the 
boundary correlation functions obtained in this paper. However, since this means we need to sum over 
intermediate spin states, the expression gets complicated. Simplifying the expression of the correlation 
functions "off the boundary" is an important problem to be considered in the future. 

Another interesting problem is to extend the analysis to other models or boundary conditions such 
as higher rank models, Felderhof model, reflecting end, etc. 
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Appendix 



We prove by induction. We show the expression holds for lattice size N and e± = — from the 
recursive relation (pTTj) . One can similarly show for e% = + from (f52"]l . Suppose (|55)l holds for lattice 
size N — 1. For ai € SL' r , one has 

^: 1 1 ' e "- ,£ ' ) ({A}\A aiJ M\u 1 ) 

= 1 A n£lj+it%^fc) y> 

det m ({\}\\ ai , m\^) ii n ^ o(A/Si ^ n ^ r+i 6(A/?) „.) ^ r ' • * J^ r 

x r_l)Ei< j<fc < s x(afc,"i) + Efc = 2( Q fc- 1 ~ r (^ + 1 )/2) + Efc =2 (^ + 1 )( Ar - fe )/2 
s 

x n ) m£j ( a °i ' a ^ ' n Eeje " ( x °i ' Aq <= ' v ^ v *) 

j=2 2<j<k<s 

xdetM({A}\{A Ql ,-- - ,Aa.},{i/}\{n,..- (A.l) 

where 

+ m \ ^ a ( A " J ' I/ i) x \ KAgjX) , A oA 

m^(A Ql , A Qj ,i/i) = -tt — r, m (A 0l , A Q . , = -r- — r. (A.2) 

We also have the follwing recursive relation [SJ for the parititon function 

Z N - 1 ({\}\\ a ,{v}\v 1 ) (-I)"" 1 A d(Ag,A a ) jM 

■?iv({A}, {i/}) a(A ct ,^i)6(A Q ,^i) Al a(A /3 ,^i)6(A /3 ,i/i) ^ a(A a , v k )b(\ a , v k ) 

detA/({A}\A a ,MVi) 

detM({A},{i/}) ' 1 J 

Combining ([5T|). (jA.lj) and (IA.3I) . one has 

j#-' e *-^({A},M) 



det M({A}, {^) 11 n - 1 a(A/J , Vj ) ft" 6(A 0; ^) 



x 



(_l)E 1 < J < fc <,x(afc.«i)+EJ= 1 (aib-l-r(6,,+l)/2)+i;j =1 ( 6 *+l)(iV-fc)/a 



X IJ II ^ efc (A Qj ,A Qfc ,^,^)detAf({A}\{A Ql ,-.- ,A as },M\K'" ,v.}) 

j=2 2<j<k<s 

c TTp=i /(Aoi,A/3) 1X0=1 rf(A,g,A ai ) s 



TTa(A Ql ,^)m^(A Ql , A Q .,!/i). (A.4) 



Since 



llfeLlKAax^fe) 



- # r (A Q1 ), 



a(A ai ,^)m + (A ai , A^-.z/i) = £ + (A Ql , A Qj , 1/1, Vj), 
a(A ttl ,^)™ _ (A Ql ,A aj ,^i) = £ (A ai , Xa^fuVj), (A.5) 



one can see that (IA.4I) is exactly the expression (1551) for ei 
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